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1.  Introduce  ion 

The  objectives  of  this  paper  are  to  present  and  analyze  a  discrete 
facility  location-allocation  model  that  is  more  general  than  similar 
models  previously  considered.  Although  the  model  is  a  deterministic  and 
static  one,  it  is  otherwise  quite  general  and  so  offers  modeling  possi¬ 
bilities  for  problem  situations  that  do  not  appear  to  have  been  formally 
treated  heretofore.  The  spirit  and  approach  of  the  paper  are  the  same  as 
those  of  an  earlier  one  [Ross  and  Soland  (1980)],  but  the  facility 
location-allocation  model  considered  herein  is  more  general  and  more 
versatile . 

To  deal  with  similarities  first,  we  feel  that  most  locational 
decision  problems,  in  both  the  public  and  private  sectors,  are  really 
multiple  criteria  problems,  and  ought  to  be  treated  as  such.  We  also 
feel  that  an  interactive  approach  to  multiple  criteria  decision  making 
(MCDM)  is  in  most  cases  the  most  appropriate  one  because  it  gets  the 
decision  maker  (DM)  involved  with  the  criteria  themselves  and  the 
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tradeoffs  possible  among  them  and  forces  him  to  directly  insert  his  pref¬ 
erences  into  the  decision  procedure.  In  our  earlier  paper  we  made  these 
points  in  connection  with  a  fairly  general  discrete  location-allocation 
model  which  we  analyzed  in  detail  with  respect  to  multiple  criteria  con¬ 
siderations.  We  discussed  at  some  length  the  various  criteria  of  cost, 
service,  and  profit  that  are  appropriate  in  public  and  private  sector 
models,  we  presented  an  introductory  look  at  MCDH  and  efficient  solutions 
and  we  outlined  a  finite  interactive  solution  algorithm  for  the  discrete 
model  that  generates  only  efficient  solutions  and  finds  them  in  a  compu¬ 
tationally  efficient  manner.  A  heuristic  procedure  based  upon  this  algo¬ 
rithm  was  subsequently  incorporated  into  a  computer  system  utilizing  a 
conversational  command  language  designed  for  flexible  user  control  [see 
Kultz,  et  al.  (1980)]. 

The  discrete  location-allocation  model  of  our  previous  paper,  like 
many  other  such  models,  does  not  directly  account  for  interactive  effects 
among  facilities  or  for  the  possibility  that  the  facilities  in  the  system 
may  be  hierarchically  structured  in  different  ways.  But  there  are 
service/distribution  systems,  such  as  multi-echelon  inventory  systems 
[Cross,  Pinkus,  and  Soland  (1981)]  and  repairable  item  support  systems 
[Cross  and  Vinkus  (1979)]  in  which  the  choice  of  hierarchical  structures 
used  is  of  great  Importance  and  partially  determines  the  facilities  to  be 
included.  And  there  are  service  systems  (e.g,,  health  care  systems)  in 
which  the  .joint  effect  of  the  facilities  utilized  determines  such  mea¬ 
sures  as  average  waiting  time  and  average  time  in  the  system.  In  general 
both  the  hierarchical  structure  and  the  interactive  effects  of  the  facil¬ 
ities  included  in  it  together  determine  some  of  the  important  criterion 
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measures  for  decision  making.  The  location-allocation  model  to  be  pre¬ 
sented  below  does  allow  for  such  generality  as  well  as  for  the  consider¬ 
ation  of  multiple  criteria,  and  it  includes  the  model  of  Ross  and  Soland 
(1980)  as  a  special  case.  It  is  based  on  the  use  of  configurations,  a 
term  we  shall  now  define.  A  configuration  is  a  specific  subset  of 

the  set  of  facilities  being  considered,  along  with  a  specific 
strategy  for  the  use  of  those  facilities.  Thus  a  configuration  has  a 
particular  hierarchical  structure  and  set  of  operating  rules,  so  two 
configurations  might  conceivably  include  the  same  facilities  and  differ 
only  in  their  structural  arrangement  or  operating  rules.  A  multi- 
echelon  inventory  system  offers  perhaps  the  best  example  of  a  configura¬ 
tion;  the  echelon  structure  and  specific  facilities  at  different  levels 
in  the  structure,  along  with  the  inventory  policies  to  be  followed  at  the 
various  facilities,  determine  the  configuration.  At  the  other  extreme, 
a  configuration  may  consist  of  a  single  facility  along  with  a  fairly 
ohvlous  and  straightforward  rationale  for  using  the  facility  to  distrib¬ 
ute  the  goods  and/or  services  for  which  it  is  responsible. 

The  other  concept  we  use  here  is  that  of  an  activity.  An  activity 
is  an  entity  that  is  required  to  be  "serviced"  by  some  particular  config¬ 
uration  in  order  for  the  system  to  function  as  intended.  An  activity  is 
to  be  carried  out  at  one  or  several  facility  locations  according  to  the 
configuration  to  which  it  is  assigned.  Thus  each  activity  is  to  be  as¬ 
signed  to  a  unique  configuration,  and  this  assignment  problem  is  what  we 
call  the  design  of  a  multiactivity  nultifacility  system.  Several  examples 
will  help  to  clarify  the  concept  of  an  activity  and  the  nature  of  the 
design  problem.  The  activities  may  be  different  products  and  the 
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configurations  multi-echelon  inventory  systems.  Since  different  products 
need  not  be  distributed  according  to  the  same  multi-echelon  structure, 
the  problem  [see  Gross,  Pinkos,  and  Soland  (1981)]  is  to  choose  an  appro¬ 
priate  configuration  for  each  product.  A  similar  problem  arises  in 
setting  up  a  repairable  item  support  system  [see  Gross  and  Pinkus  (1979)]; 
<acli  different  repairable  component  constitutes  an  activity  and  each  con¬ 
figuration  is  a  subset  of  facilities  at  which  to  perform  the  repairs, 
along  with  rules  to  apportion  the  population  of  that  component  among  the 
facilities  in  the  given  subset.  An  activity  may  be  a  particular  type  of 
health  care  or  educational  unit  (such  as  a  cardiac  unit)  or  a  par¬ 
ticular  type  of  emergency  service  vehicle  (such  as  an  ambulance  or  small 
fire  engine)  and  the  configurations  are  then  different  subsets  of  facil¬ 
ities  at  which  to  locate  these  units  or  vehicles. 

Having  indicated  the  similarities  and  differences  between  the 
present  paper  and  our  previous  one,  we  remark  that  the  reader  is  not 
expected  to  be  familiar  with  the  previous  results.  The  presentation  here 
is  self-contained,  but  some  aspects  of  it  will  be  covered  briefly  in 
order  to  avoid  unnecessary  repetition. 

In  the  remainder  of  this  section  we  briefly  discuss  the  different 
criteria  which  ought  to  be  considered  by  those  responsible  for  the  design 
of  raultiactivity  multifacility  systems.  In  Section  2  we  formulate  a  dis¬ 
crete  choice  model  for  the  design  of  multiactivity  nultifacility  systems 
and  we  quantify  the  various  criteria  previously  discussed  in  terms  of  the 
variables  and  parameters  of  this  model.  We  also  show  that  our  previous 
model,  and  hence  a  wide  class  of  location-allocation  models,  is  obtained 
as  a  special  case  of  the  present  model.  In  Section  3  we  briefly  discuss 
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the  main  approaches  to  MCDtt,  present  a  general  characterisation  of  the 
efficient  (or  nondominated)  solutions  of  a  MCDM  problem,  and  specialize  it 
to  the  discrete  model  of  Section  2  for  the  design  of  multiactivity  multi- 
facility  systems.  In  Section  4  we  introduce  a  0-3  linear  optimization 
model  called  the  team  generalized  assignment  problem  (T-GAP)  and  show  that 
parametric  solution  of  a  T-GAP  will  yield  all  efficient  solutions  of  the 
multiactivity  multifacility  design  problem  with  multiple  criteria.  We 
also  discuss  there  the  numerical  solution  of  T-GAPs.  Rather  than  attempt¬ 
ing  to  find  all  efficient  solutions  of  the  design  problem,  however,  we 
advocate  an  interactive  approach  and  briefly  describe  in  Section  5  an 
interactive  branch-and-bound  algorithm  that  solves  a  finite  sequence  of 
T-GAPs  in  order  to  determine  a  solution  to  the  design  problem  that  is 
most  preferable  to  the  DM.  Some  final  remarks  conclude  the  paper. 

Now  we  turn  to  a  discussion  of  the  criteria  to  be  considered. 

As  indicated  above,  the  criteria  naturally  fall  into  the  categories  of 
cost,  service,  and  profit,  but  there  exist  several  distinct  criteria  in 
each  of  these  categories.  Most  of  them  have  been  discussed  by  Ross  and 
Soland  (1980,  Sections  1  and  5),  so  we  shall  merely  cite  those  here. 

With  respect  to  costs,  we  deal  with  three  classes:  investment  cost, 
operating  cost,  and  discounted  cost.  Within  each  of  these  classes  we 
differentiate  between  the  fixed  cost ,  the  incremental  (or  variable)  cost, 
and  the  total  east,  which  is  the  sum  of  the  fixed  and  incremental  com¬ 
ponents.  The  fixed  cost,  in  each  of  the  three  classes,  is  assumed  to  be 
independent  of  the  size  of  the  facility  concerned  and  the  services  of¬ 
fered  there.  The  incremental  cost  at  a  particular  facility,  on  the  other 
hand,  is  completely  dependent  on  the  facility  size  and  services  offered. 
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But  in  the  context  of  multiactivity  raultifacility  systems  there  are  also 
incremental  costs  that  need  not  necessarily  be  allocated  to  a  particular 
facility.  For  example,  in  a  multi-echelon  inventory  system  the  total 
annua]  inventory  cost  for  a  particular  product  may  be  of  more  interest 
than  the  various  inventory  costs  at  the  respective  facilities.  Indeed, 
the  inventory  policy  followed  by  the  system  (for  this  particular  product) 
may  be  chosen  to  minimize  this  quantity.  Thus,  while  fixed  costs  will 
always  be  attributed  to  specific  facilities,  it  will  be  convenient  to 
deal  with  incremental  costs  that  are  attributed  to  configurations  as 
well  as  those  that  are  attributed  to  specific  facilities. 

We  distinguish  between  investment  and  operating  costs  on  the  basis 
of  when  they  occur — before  or  after  the  facilities  are  in  operation. 

Total  discounted  cost  is  in  appropriate  weighted  sum  of  investment  and 
operating  costs. 

Now  we  turn  to  service  criteria.  These  may  differ  significantly, 
depending  on  the  context  of  the  problem.  For  a  system  providing  services 
to  the  public,  the  total  demand  served  and  the  demand  which  can  be  served 
within  a  specified  time  (or  distance)  are  important  measures.  Also  im¬ 
portant  is  the  average  response  time  or  average  distance  traveled.  The 
average  time  spent  waiting  for  service  and  average  total  trip  time  are 
other  important  ones.  For  a  multi-echelon  inventory  system  the  overall 
availability  and  mean  turnaround  time  of  orders  may  be  the  most  appropri¬ 
ate  service  criteria,  and  similar  ones  may  be  used  for  repairable  item 
support  systems. 

Criteria  associated  with  profit  are  important  in  the  private  sec¬ 
tor;  both  annual  profit  and  total  discounted  profit  are  appropriate  ones. 
Total  annual  revenue  may  also  be  of  interest,  and  we  include  it  as 
another  criterion  associated  with  profit. 
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Although  certain  criteria  are  to  be  ninitaized  while  others  are  to 
be  maximized,  it  will  hereafter  be  convenient  to  talk  of  minimizing  all 
criteria.  It  will  thus  be  necessary  to  use  Lhe  negatives  of  certain 
criteria;  e.g.,  the  negative  of  demand  served  and  the  negative  of  profit. 

•  •  A  Discrete  Location  Model 

The  nultiactivity  multifacility  design  problem  is  to  assign  each 
activity  to  a  specific  configuration  in  a  manner  that  respects  capacity 
constraints  at  the  facilities  utilized;  for  obvious  reasons  we  do  not 
specify  an  objective  at  this  point.  We  suppose  there  are  N  distinct 
activities  to  be  assigned  and  I.  distinct  sites,  at  some  subset  of 
which  facilities  will  be  (or  already  are)  established.  We  number  from 
1  to  M,  the  configurations  to  which  activity  j  nay  be  assigned  and 
note  that  configuration  i  for  activity  j  need  not  bear  any  relation 
to  configuration  i  for  activity  j'  .  We  assume  that  configurations 
1  to  are  specified  for  each  j  and  do  not  deal  here 

with  the  question  of  generating  the  sets  of  configurations  to  be  consid¬ 
ered.  This  is  not  a  trivial  problem,  and  we  shall  have  some  more  to 
say  about  it  in  Che  final  section. 

To  simplify  our  notation  and  terminology  we  let  M  *  Max  M  and 

j  3 

suppose  tnat  each  of  the  N  activities  is  to  be  assigned  to  one  of  M 
conf igerut ions ;  as  will  be  seen,  it  is  a  simple  matter  to  prevent  as¬ 
signment  of  activity  j  to  configuration  i  if  i  exceeds  . 

For  h«l,...,t,  ,  i-l,,..,M  ,  and  J-1,...,N  ,  let 
-  1  if  a  facility  is  established  at  site  b, 
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x^  *  1  if  activity  j  is  assigned  to  configuration  i. 


0  otherwise; 


d^„  =  the  usage  of  the  facility  at  site  h  due  to  the  assignment 


of  activity  j  to  configuration  i; 


a,  ,b,  *  the  minimum  and  maximum,  respectively,  total  usage  per¬ 
il  n 


mitted  at  site  h. 


Feasible  solutions  of  the  multiactivity  multifacility  design 


problem  are  then  those  choices  of  the  z.  and  x.,  which  satisfy  the 


following  constraints: 


N  H 


<  Ji  ‘‘hij-lj  <  Vh  •  f°ral‘  '"l . L  ■ 


l  x,.  »  1  ,  for  all  j=l,...,N  , 


Vx  =  0  or  1  ,  for  all  h-l,...,L;  i«l,...,M;  j-l,...,N  .  (3) 


Clearly,  0  <  a,  <  b,  for  all  h  ,  but  it  is  quite  permissible  to  have 
n  n 


a.  •*=  0  and/or  b.  °  B  (a  very  large  number  which  represents  +*>)  to  in- 
h  n 


dicate  the  lack  of  limitations  on  total  usage  at  facility  h  .  Each 


dhlj  is  nonnegative,  but  clearly  d^  ■  0  if  configuration  i  does 


not  involve  the  use  of  facility  h  for  activity  j  .  Clearly,  the 


d)itj  must  be  so  defined  that  the  quantity  ^  ^  dhijXij  is  3  meanin8~ 


ful  measure  of  the  usage  of  facility  h  .  Determination  of  d^^  (and 


of  other  parameters  to  be  introduced)  is  not  necessarily  simple,  and 


may  indeed  involve  the  solution  of  an  optimization  problem.  For  example, 


In  the  case  of  multi-echelon  inventory  systems  [Gross,  Pinkus,  and 


Soiand  (1981)1,  configuration  i  represents  a  particular  hierarchical 


structure  along  with  the  minimum  cost  inventory  policy  (for  a  specified 


■  . . . 
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product)  using  that  structure.  To  prevent  the  assignment  of  activity  j 
to  configuration  i  if  i  >  we  may  set  dj^  =  +°°  for  all  h  . 

Now  we  turn  to  the  criterion  functions  to  be  used  with  this  model. 
We  use  the  index  k  to  distinguish  various  criteria.  For  investment, 
operating  or  discounted  cost,  the  quantity 


^  1 hkZh 

h=l 


(4) 


is  the  fixed  cost,  whereas  we  take 


N  M 


l  1  C  x, 
j-i  i-1  iJU  ij 


(5) 


and 


Jl  ''**»  +  ih  X  (6> 

as  the  incremental  cost  and  total  cost,  respectively.  The  and 

are  nonnegative  constants  whose  determination,  as  noted  before, 
nay  not  be  simple. 

For  service  and  profit  criteria  the  exact  forms  of  the  criterion 
functions  depend  on  the  context  as  well  as  on  the  criterion,  and  the 
reader  is  referred  to  Ross  and  Soland  (1980,  Sections  2  and  5)  for  de¬ 
tailed  treatment  of  a  number  of  specific  cases.  What  is  important,  how¬ 
ever,  is  that  almost  all  reasonable  criterion  functions  can  be  written 
in  the  form 


L  N  M 

£  fhWZh  +  1  l  S-MkX-H  * 

h“l 


hk  h  ,L.  “ijk"ij 

j=l  i=l  J 


(7) 


where  the  f ^  and  are  appropriately  defined  constants.  In  gen¬ 

eral  we  wish  to  minimize  each  of  several  expressions  of  the  form  (7), 
subject  to  the  constraints  (1)  -  (3). 
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The  multiactivity  multifaciltty  design  problem  has  not  previously 
been  considered  in  the  generality  presented  here,  especially  with  re¬ 
spect  to  multiple  criteria,  but  its  history  goes  back  at  least  to  the 
dissertation  of  Pinkus  (1971)  and  the  published  version  of  Pinkus  (1975). 
Pinkus,  Gross,  and  Soland  (1973)  considered  the  problem  in  more  general¬ 
ity  than  the  inventory  context  and  concentrated  on  a  specialized  branch- 
and-hound  solution  scheme.  All  three  of  these  papers  dealt  with  the 
uncapacitated  case  (all  *  0  and  all  =  B)  ,  as  did  the  applica¬ 

tion  of  Gross  and  Pinkus  (1979).  The  capacitated  case  was  used  us  a 
model  by  Gross,  Pinkus,  and  Soland  (1981),  and  Chhabra  and  Soland  (1980) 
and  Chhabra  (1981)  have  presented  a  non-LP  branch-and-bound  algorithm 
and  computer  code  for  this  case.  Ail  the  papers  cited  here  dealt  with  a 
single  criterion  function  of  the  form  (7). 

We  indicated  earlier  that  the  present  model  includes  that  of  Ross 
and  Soland  (1980)  as  a  special  case.  To  see  this,  let  each  activity 
represent  a  client  group  and  let  each  configuration  represent  a  single 
facility  site  to  which  client  groups  may  be  assigned.  Thus,  *  M  *  L 
for  all  ;  we  keep  the  natural  correspondence  between  configu¬ 

rations  and  sites  so  that  the  index  h  is  replaced  by  i  and  L  is 
replaced  by  M  .  Also,  is  zero  unless  h=i  ,  in  which  case  we 

call  it  d^  .  Hence  (1)  reduces  to 


atzi  <  i  d 


j-i 


ijxij  <  Vi  * 


for  all  i=l , . . . ,M  , 


(8) 


which  is  Equation  (1)  of  Ross  and  Soland  (1980),  It  is  also  easy  to 
see  that  each  criterion  function  of  the  general  form  (7)  reduces  to  the 
form  used  in  the  earlier  model. 
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3.  Multiple  Criteria  Optimization  and 
Efficient  Solutions 

Most  techniques  for  deterministic  MCDM  involve  some  combination 
of  these  three  approaches:  construction  and  use  of  a  vulue  function, 
determination  of  the  set  of  efficient  solutions,  and  use  of  an  inter¬ 
active  algorithm  [see  Soland  (1979)  for  an  overview  of  these  approaches 
and  a  number  of  references] .  In  a  general  treatment  of  MCDM  we  take  X 
as  the  set  of  feasible  solutions  and  assume  there  are  p  criterion 
functions,  ,  all  to  be  minimized.  We  write  f(x)  for  the 

vector  (f^(x),. . . ,fp(x))  ,  and  phrase  the  MCDM  problem  as 


'Minimize”  f(x) 

(9a) 

s.t.  x  e  X  , 

(9b) 

Let  V  =  {y  c  RP  |  y  =  f  (x)  for  some  xeX)  be  the  set  of  feasible  cri¬ 
terion  vectors. 

The  value  function  approach  is  based  on  the  result  that,  under 

certain  hypotheses  concerning  the  preference  structure  of  the  DM,  there 

exists  an  isotone  decreasing  value  function  v  defined  on  a  set  con- 

1  2 

talning  Y  such  that  the  DM  prefers  y  e  Y  to  y  e Y  if  and  only  if 
1  2 

v(y  )  >  v(y  )  .  If  the  DM's  value  function  v  is  given,  a  most  desir¬ 
able  solution  to  his  MCDM  problem  is  clearly  one  that  is  optimal  for 
problem  (10) : 

Maximize  v(f(x))  (10a) 

s.t.  x  e  X  .  (10b) 

The  major  difficulty  in  practice  is  the  determination  of  v  ;  see 
Keeney  and  Raiffa  (2976)  for  an  extensive  treatment. 
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At  the  opposite  extreme  from  the  use  of  a  value  function  is  the 
determination  of  the  set  of  efficient  solutions  of  (9);  this  demands 
practically  no  preference  specification  on  the  part  of  the  DM.  We 
start  with  the  definition  of  a  dominated  solution: 


Definition:  t  X  is  a  dominated  solution  of  (9)  if  there  exists 

x  t:  X  such  that  f(x)  <  f(x^)  but  f(x)  f  f(x^*)  . 


An  Ax  that  is  not  dominated  is  said  to  be.  nondominated  or 
efficient;  we  let  X£  denote  the  set  of  efficient  solutions  of  (9). 

A  second  approach  to  the  resolution  of  (9)  is  to  present  to  the  DM  the 
efficient  set  Y  =  f(X„)  «  the  set  of  efficient  solutions  in  criterion 
space.  A  result  which,  at  least  in  theory,  makes  it  easier  to  determine 
X„  and  Y„  Is  based  on  the  auxiliary  problem  (P,  ) : 

K  h  A  f  (X 


P 

Minimize  Xf (x)  *  £  X  f  (x) 

k-1 


s.t.  x  e  X  , 

f(x)  <  a  , 


where  X  >  0  and  a  e  Rp 


Then  we  have 


) 


Lemma  1  [Soland  (1979)]:  For  arbitrary  X >  0  ,  x° c  if  and  only  if 
x°  is  optimal  in  (P,  )  for  some  a  £  RP  . 

A 

Thus  parametric  solution  of  problem  (P,  )  will  yield  all  elements  of 

A  9  Ot 

X  and  '  .  Note  that  this  result  is  independent  of  the  nature  of 

tli 

the  constraint  set  X  and  the  criterion  functions  f  and  so  applies, 
e.g.,  to  discrete  problems.  It  will  be  used  below  to  characterize  the 
efficient  solutions  of  our  multiactivity  nultifacility  design  model 


-  12  - 


T-444 


(1)  -  (2)  a«d  (7).  But  first  we  introduce  the  interactive  approach  to 
MCDM,  with  which  we  shall  deal  exclusively  in  Section  5. 

In  the  interactive  approach  the  DM  works  directly  with  an  ana¬ 
lyst  and/or  an  interactive  computer  program.  The  DM  supplies  informa¬ 
tion  which  helps  to  reveal  his  preferences  while  the  analyst  and/or 
computer  program  is  responsible  for  generating  solutions  for  the  DM  to 
consider  and  evaluate.  An  interactive  procedure  should,  after  a  rea¬ 
sonable  amount  of  effort  on  the  part  of  the  DM,  yield  a  solution  with 
which  he  is  "satisfied.11  It  should  be  an  efficient  solution,  and  it  is 
desirable  that  it  be  one  that  maximizes  the  DM's  (unknown)  value  func¬ 
tion.  Most  interactive  procedures  make  certain  convexity  or  linearity 
assumptions  about  the  criterion  functions  and  constraint  set  and/or 
the  DM's  value  function.  In  Section  5  we  describe  an  interactive  algo¬ 
rithm  that  does  not  have  these  limitations  and  which  we  hope  to  be 
computationally  efficient  for  our  design  problem. 

Now  we  return  to  the  multiactivity  multifacility  design  problem 
of  Sectio  2  and  its  set  XE  of  efficient  solutions.  We  write  (x,z) 
for  the  vector  of  x^  and  z^  values  that  identify  a  feasible  solu¬ 
tion  of  (1)  -  (3).  We  suppose  that  among  the  many  possible  criterion 
functions  discussed  exactly  p  ,  all  having  the  algebraic  form  (7),  are 
selected  by  the  DM  as  relevant  for  the  design  decision.  Then,  by  vir¬ 
tue  of  I.enma  1,  we  have  the  following  result: 

Corollary  1.1;  for  arbitrary  A>0  ,  (x,z)  e  Xp  if  and  only  if  (x,z) 
Is  optimal  in  the  following  auxiliary  problem  (P,  )  for  some  a  v.  RP  : 
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P  L  N'  M 


s.  t  , 


N  M 


<U) 


(PX.a> 


Vh  «  J.  J,  dhijxii  *  Vh  *  h_1 . L  1  (1) 

3=1  i=l  J 


I  x-  •  =  1  » 

l  11  ’ 


i=l 


L  N  M 

^  fhk2h+  2  2  S'-’X- 


h»l 


J-Iiii  * 


,  (2) 


a  ,  k=l,...,p  ,  (12) 


Zh’Xij  “  0  or  X;  h=l,...,L;  1=1,..., M;  j=l,...,N  ,  (3) 

where  >  0  ,  c  R1  ,  and  the  and  are  appropriately 

chosen  constants.  Note  that  the  objective  function  (11)  of  (P^  a)  can 
be  written  as 

L  N  M 

I  fh(A)zh  +  $  l  sii(A)xij  * 

h-1  n  n  j-1  i-1 


dl') 


where  f.  (X)  T,  X.  f, ,  and  s.,(A)  »  J.A.  s...  .  In  the  next  section 

h  k  hk  lj  Lk  k  ijk 

we  shall  relate  problem  (P,  )  to  a  discrete  optimization  model  we  call 

A  j  CX 

the  team  generalized  assignment  problem  (T-GAP). 


1 .  The  Team  Generalized  Assignment 
Prob lem  (T-GAP) 

The  generalized  assignment  problem  (GAP)  ,  sot  and  Soland 

(1975)  and  Ross,  Soland,  and  Zoltners  (1980)]  is  a  discrete  optimiza¬ 
tion  model  whjch  nay  be  interpreted  as  a  problem  of  assigning  tasks  to 
agents  at  ninir  ira  cost,  subject  to  upper  and  lower  limits  on  the  amount 
of  a  resource  expended  by  each  agent.  We  here  generalize  that  model  so 
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that  tasks  ate  assigned  to  teams,  each  team  being  composed  of  one  or 
more  agents  and  having  a  defined  strategy  for  how  the  agents  work  to¬ 
gether  to  perform  a  given  task.  A  mathematical  formulation  of  the 
T-GAP  is: 


Minimize 

l  l  rOijXij 

jeJ  id  3  13 

(13) 

3.  t. 

a,  <  y  y  r.  ,  .x  , 

h  ~  j£j  id  hij  lj 

Z  bh  * 

for  all  h  C  H  , 

(14) 

l  xi-j  =  1  * 
id  J 

for  all  j  f.  J  , 

(15) 

x. .  =  0  or  1  , 
ij 

for 

all  i  e  I,  j  e  J  . 

(lb) 

H  E  {1,2, 

is  the  set  of 

agent 

indices,  1  E  {1,2,. 

•  •  j  III  } 

is  the  set  of  team  indices,  J  =  {l,2,...,n}  is  the  set  of  task  in¬ 
dices,  ro£j  is  t^ie  cost  incurred  if  task  j  is  assigned  to  team  i  , 
r^j  >  0  is  the  amount  of  a  resource  required  by  agent  h  to  help 
perform  task  j  as  a  member  of  team  i  ,  and  >  0  and  b^  >  0 
are,  respectively,  the  minimum  and  maximum  amounts  of  the  resource 
that  nay  be  expended  by  agent  b  .  The  decision  variables  are  inter¬ 
preted  as 

x^  =1  if  task  j  is  assigned  to  team  i, 

*  0  otherwise. 

We  note  that  rhij  =  ®  if  agent  h  is  not  a  member  of  team  i  ,  and 
we  remark  that  one  need  not  actually  consider  the  same  number  or  compo¬ 
sition  of  teams  for  all  the  tasks  (so  a  team  is,  intentionally,  similar 
to  a  configuration).  This  generalization  is  dealt  with  as  in  Section  2 
We  now  proceed  to  show  that  problem  (P,  )  can  be  written  as  a 

A  jOt 

T-GAP  of  size  H  -  L+p,  m  ■  M  ,  and  n  =  N+L  ,  In  the  T-GAP  formula¬ 
tion  the  x.,  are  defined  as  in  (P.  )  for  1=1,... ,M  and  j=l,...,N 
lj  A  ,a 
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The  variable  2^  af  (P^  a)  is  represented  by  ,  so  we  define 

x,  v,.,  =1  if  a  facility  is  established  at  site  h, 

J  ,N+n 


and 


=  0  otherwise: 


x„  *0  if  a  facility  is  established  at  site  h, 

l  ,N+h 


=  1  otherwise. 

Thus  x^  =  i  for  h=l .  All  other  for  j  >  N 

must  be  zero.  We  now  give  a  full  stipulation  of  the  T-GAP  parameters. 

For  h  <  L  the  T-GAP  parameters  a,  and  b,  are  defined  as  in  (P.  ). 
-  h  h  A, a 

For  h  -  L-Pk  ,  k=l,...,p  ,  we  distinguish  two  cases  which  correspond  to 

whether  the  k th  natural  criterion  function  is  one  to  be  minimized  (such 

as  cost)  or  maximized  (such  as  demand  served).  In  the  former  case 

>  0  and  Sijk  ®  f°r  h,t,j  and  ai^  -■  0  also.  Then  we  set 

a.  =  0  and  b.  =  a  for  h  =  L+k  .  In  the  latter  case  f,_  <  0  and 
h  h  k  hk  • 

<  0  for  all  h,i,j  and  <  0  also.  Then  for  (12)  we  have 
the  equivalent  constraint 


L 

N 

KJ  -  J.  |fhkizh 

h=l 

+  ^ 
.1-1 

so  we  set 

a,  *  |u,  I  and  b  =  B 
h  1  k'  h 

for  h 

roij 

= 

slj(A)  , 

i=l 

r0,l,N+h 

= 

fh(A)  r0,2,N+h  =  0 

,  h-1 

rQ, i,N+h 

= 

00  ,  for  1  >  2  and 

h=l 

Chij 

s 

dhij 

h-1 

rh,2,N+h 

s 

bn  * 

h-1 

rhlj 

GZ 

0  otherwise, 

h=l 

(17) 
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L+k,i,j  ‘Sijk‘  ’ 

i=l, . . . ,M; 

j=l,. 

..,N; 

=  If  1 

L+k ,  1 ,  N+h  1  hk1  ’ 

h“ 1 , - • • , L; 

k=l, . 

•  •  *P 

,  „  =  o  , 

h=l,. . . ,L; 

k=l , . 

. .  ,p 

L+k, 2, N+h  ’ 

•  ,  =  00  otherwise, 

bij 

Thus  problem  (P,  )  is  equivalent  to  a  T-GAP,  and  hence  by 

A  y  Ot 

Corollary  1.1  all  efficient  solutions  of  the  multiactivity  multifacil¬ 
ity  design  problem  with  multiple  criteria  can  be  found  by  parametric 
solution  of  this  equivalent  T-GAP.  Rather  than  attempting  to  find  all 
efficient  solutions,  however,  we  advocate  an  interactive  approach  that 
utilizes  the  equivalent  T-GAP  for  computational  efficiency.  We  will 
return  to  this  point  in  the  next  section,  but  first  we  address  the 
question  of  efficient  numerical  solution  of  the  T-GAP. 

The  T-GAP  is  a  somewhat  specialized  version  of  the  multiple 
choice  integer  program  (MCIP) ,  the  specialization  being  that  the  con¬ 
straint  coefficients  in  (14)  are  assumed  to  be  nonnegative.  The 

constraints  (15)  and  (16)  together  are  termed  multiple  choice  con¬ 
straints.  Bean  (1980)  discusses  the  MCIP  at  length,  provides  refer¬ 
ences,  and  presents  a  new  algorithm  for  its  solution.  He  cites  encour¬ 
aging  computational  results  for  some  problems  with  up  to  400  variables. 
Another  approach  for  solving  the  MCIP  is  given  by  Clover  and  Mulvey 
(1979). 

It  seems  to  us  that  a  combination  of  approaches  for  the  MCIP 
will  eventually  yield  the  most  efficient  algorithm  for  the  T-GAP. 
Besides  the  additive  approach  of  Bean  (1980),  some  combination  of 
Lagrangian  relaxation  and  surrogate  constraints  [see  Karwan  and  Rardin 
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(1979)  and  Dyer  (I960) J  will  be  important  in  generating  tight  lower 
bounds  for  a  branch-and-bound  algorithm.  Lagrangian  relaxation  has  al¬ 
ready  proved  very  successful  in  solving  certain  special  classes  of 
MCIPs  that  are  location  problems  [see  Fisher  (1981)  for  a  discussion}; 
both  subgradient  optimization  [see  Crowder  (1976)]  and  multiplier  ad¬ 
justment  methods  [see  Fisher  (1981)]  have  proved  successful  in  provid¬ 
ing  the  tight  Lagrangian  bounds  desired. 

Chhabra  (198i)  and  Chhabra  and  Soland  (1980)  have  directly 
attacked  the  multiactivity  multifacility  design  problem  with  a  single 
criterion  function.  They  treat  the  case  with  all  =  0  and  use  a 
combination  of  Lagrangian  relaxation  and  the  additive  approach. 

5.  Interactive  Solution  of  the  Multiple 
Criteria  Location  Model 

In  Ross  and  Soland  (1980,  Section  4)  we  presented  an  interactive 
satisficing  algorithm  for  MCDM  that  terminates  after  generating  a  fi¬ 
nite  number  of  efficient  solutions.  It  does  not  generate  any  dominated 
solutions.  In  it  the  DM  is  asked  to  make  some  binary  comparisons  of 
vectors  in  the  efficient  set  Yg  .  He  is  also  asked  to  vary  the 
weights  and  to  gradually  decrease  the  "satisfaction  levels" 

that  appear  in  problem  (P,  ) .  That  algorithm  may  of  course  be  applied 

A  f  (X 

in  the  current  context  without  alteration,  and  in  that  case  each  of 
the  finite  number  of  optimization  problems  to  be  solved  is  a  T-CAP; 
only  A  and  a  change  from  one  problem  to  the  next. 

Hultz,  ot  at.  (1980)  subsequently  incorporated  a  heuristic  pro¬ 
cedure  based  upon  the  above  algorithm  into  an  interactive  computer 
system  utilizing  a  conversational  command  language.  That  same  computer 
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system,  with  rainor  alteration,  could  be  modified  to  deal  with  the 
present  multiple  criteria  location  model,  and  would  then  he  required 
to  solve  a  T-CAF  at  each  Iteration. 

lhe  two  interactive  algorithms  just  mentioned  nav  be  criticized 
as  being  too  demanding  of  the  DM  in  that  they  require  him  to  provide  a 
,-equence  of  X  and  a  vectors.  Equally  important,  it  is  not  clear 
how  the  X  and  rx  vectors  can  be  related  to  the  DM's  preference 
structure,  and  so  the  algorithms  do  not  necessarily  yield  a  decision 
choice  that  is  in  some  sense  "optimal"  for  the  DM.  In  contrast,  the 
interactive  algorithm  of  Marcotte  and  Soland  (1981),  to  be  briefly 
described  here,  only  requires  the  DM  to  make  binary  comparisons  of 
vectors  in  criteria  space  and  finally  yields  a  decision  choice  that 
solves  problem  (10),  i.e.,  maximizes  the  DM's  (unknown)  value  function. 

The  algorithm  of  Marcotte  and  Soland  (1981)  uses  a  branch-and- 
bound  approach  to  solve  the  following  problem,  which  is  equivalent 
to  (10); 

Maximize  v(y)  (18a) 

s.t.  y  e  Y  .  (18b) 

The  feasible  set  Y  ,  and  subsequently  subsets  of  Y  ,  are 

separated  into  subsets  by  appending  to  the  requirement  y  c  Y  addi¬ 
tional  constraints  of  the  form  y^  <  ,  one  for  each  subset  gener¬ 

ated.  The  otk  may  be  interpreted  as  satisfaction  levels,  and  their 

values  are  obtained  from  an  efficient  point  y  e  Y_  .  We  illustrate 

h 

this  separation  process  with  the  inir  .1  subset  YQ  =  Y  .  An  efficient 
point  y°  is  found  as  the  solution  of  the  problem 
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Minimize  j  c,  y 

(1 9a) 

k3 1 

s.t.  y  *0  . 

(19b) 

where  the  are  positive  constants;  Y^  is  then  separated  into  p 

subsets  ,  k=l,...,p  ,  defined  as 

Y.  =  Y,  A  {v  |  v  <  y1)  ,  k=l,..,,p  .  (20) 

K  ij  K  K 

For  discrete  problems,  such  as  those  we  deal  with  here,  the  strict 
inequality  in  (20)  causes  no  difficulty. 

In  an  ordinary  branch-and-bound  algorithm  one  would  compute  an 
upper  bound  on  v  over  the  subset  Y^  for  each  k=l,...,p  .  All  sub¬ 
sets  whose  upper  bounds  dla  not  exceed  the  value  of  the  best  known 
feasible  solution  (the  incumbent  solution  y*  )  would  be  discarded, 
ana  branching  would  continue  from  one  of  the  remaining  subsets.  In 

this  interactive  algorithm  an  upper  bound  is  the  value  of  v  asso- 

k  1c  1c  k 

ciated  with  the  ideal  vector  6  *  (gj,...,g  )  ,  where  g^  is  the  op¬ 
timal  value  of  the  problem 

Minimize  y^  (21a) 

s.t.  y  e  Yk  ,  (21b) 

k  k 

£“l,...,p  .  Note  that,  in  general,  g  t  and  g' '  t  Y  .  At  this 

point  of  the  algorithm  the  incumbent  solution  is  y*  =  ,  and  it  is 

0  k 

necessary  to  compare  v(y  )  with  the  upper  bounds  v(g  )  ,  k=l,...,p  . 

As  v  is  unknown,  such  comparisons  are  not  possible.  Instead,  the  DM 

0  k 

is  asked  to  compare  the  two  vectors  y  and  B  and  to  stac6  which  he 

prefers;  this  comparison  is  equivalent  to  one  between  v(y^)  and 

.  -k. 

v(g  )  . 
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There  are  other  comparisons  required  by  this  brnnch-aud-baund 
algorithm  that  are  performed  by  the  PM.  He  must  compare  y*  with 
other  efficient  solutions  generated,  and  he  musL  sometimes  arrange 

the  subsets  generated  at  a  separation  according  to  decreasing 

preference  of  their  ideal  vectors  B^  .  This  latter  task  can  be  done, 

■i the  DM  so  desires,  through  a  series  of  binary  comparisons. 

Thus  the  algorithm  only  requires  the  DM  to  make  binary  compari¬ 
sons  of  vectors  in  criteria  space.  It  generates  a  finite  sequence  of 
feasible  solutions,  all  of  which  are  efficient,  and  it  terminates  after 
a  finite  number  of  separations.  Assuming  that  the  DM's  responses  arc 
consistent  with  his  preference  structure,  the  incumbent  solution  at  the 
termination  of  the  algorithm  solves  problem  (18),  i.e.,  is  a  most  pre¬ 
ferred  choice.  For  each  separation  the  algorithm  is  required  to  solve 
one  problem  of  the  form  (19)  and  (up  to)  p  problems  of  the  form  (21). 
As  we  have  seen  In  Section  4,  these  problems  may  be  solved  as  T-GAPs; 
only  the  X  and  a  vectors  change  from  one  problem  to  another.  With 
an  efficient  method  for  solving  such  T-GAPs  the  algorithm  of  Marcotte 
and  Soland  may  become  a  practical  vehicle  for  dealing  with  the  wide 
class  of  locational  decisions  that  may  be  modeled  as  multiactivity 
multifacility  design  problems  with  multiple  criteria. 

6 .  Concluding  Remarks 

The  locational  decision  model  presented  here  is  more  general 
than  other  static  deterministic  location  models  and  allows  for  the  in¬ 
corporation  of  interactive  effects  among  facilities  and  hierarchical 
structures  ol  the  facilities.  For  models  of  significant  size,  however, 

-  21  - 


T-444 


fairly  difficult  0-1  optimization  problems  result.  These  can  be 
reduced  to  the  form  of  T-GAPs,  for  which  we  hope  that  computationally 
efficient  computer  codes  will  become  available.  Work  is  definitely 
needed  on  efficient  algorithms  for  the  T-GAP. 

The  coefficients  f,  ,  s, ,  and  d,  ,  ,  of  our  model  depend 
k  ijk  *  hij 

:  oavily  upon  the  context  of  the  problem  and  may  be  difficult  to  deter¬ 
mine  in  some  cases.  Monte  Carlo  simulation  may  even  be  required  to 
estimate  them  in  certain  applications.  Clearly,  then,  good  submodels 
are  needed  to  generate  the  model  coefficients  for  specific  problems. 

It  is  also  necessary  to  generate  ’’good"  configurations  to  be 
considered  for  the  various  activities  to  be  assigned.  In  most  cases 
it  will  not  be  possible  to  consider  all  possible  configurations,  so  it 
will  partially  fall  upon  the  tnodeler/analyst  to  insure  that  the  subset 
of  configurations  explicitly  considered  by  the  model  is  sufficient  ;y 
"rich."  Some  suboptimization  may  inevitably  result  from  the  interac¬ 
tion  of  the  capacity  constraints  and  the  inability  to  include  all 
possible  configurations;  this  difficulty  is  worthy  of  examination  and 
may  he  surmountable  in  some  cases. 

We  have  stressed  here  an  interactive  approach  to  an  MCDM  prob¬ 
lem,  but  we  remark  that  the  model  and  T-GAP  equivalency  can  be  useful 
in  generating  efficient  solutions  with  guaranteed  maximal  values  of  the 
criterion  functions  even  if  no  explicit  interactive  procedure  is  used, 
"ven  if  one  criterion  is  of  paramount  importance,  the  model  can  be  used 
to  impose  constraints  on  a  number  of  less  important  criterion  functions, 
e.g.,  an  upper  limit  may  be  set  on  the  total  number  of  facilities. 
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